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We make a detailed study on the Higgs-boson decay width T(H —>• gg) up to a®-order under 
the minimal momentum space subtraction (mMOM) scheme. A major uncertainty of a finite-order 
perturbative QCD prediction is the perceived ambiguity in setting the renormalization scale. In 
the present paper, to achieve a precise pQCD prediction without renormalization scale uncertainty, 
we adopt the Principle of Maximum Conformality (PMC) to set the renormalization scale of the 
process. The PMC has solid theoretical foundation, which is based on the renormalization group 
invariance and utilizes the renormalization group equation to fix the renormalization scale of the 
process. The key point of the PMC applications is how to correctly set the {/3i}-terms of the 
process such that to achieve the correct a s -running behavior at each perturbative order. It is 
found that the ambiguities in dealing with the {/3i}~terms of the decay width T(H —► gg) under 
the MS-scheme can be avoided by using the physical mMOM-scheme. For the purpose, as the first 
time, we provide the PMC scale-setting formulas within the mMOM-scheme up to four-loop level. 
By using PMC, it is found that a more reliable pQCD prediction on V(H —► gg) can indeed be 
achieved under the mMOM-scheme. As a byproduct, the convergence of the resultant pQCD series 
has also been greatly improved due to the elimination of renormalon terms. By taking the newly 
measured Higgs mass, Mh = 125.09 ± 0.21 ± 0.11 GeV, our PMC prediction of the decay width 
is, r (H —> gg)| m MOM,PMC = 339.3 ± 1.7~^ 3 I keV, in which the first error is from the Higgs mass 
uncertainty and the second error is the residual renormalization scale dependence by varying the 
initial renormalization scale p r € [Mh/2,4Mh], 

PACS numbers: 14.80.Bn, 12.38.Bx, 12.38.Cy 

I. INTRODUCTION 

Both the ATLAS and the CMS collaborations at the 
Large Hadron Collider (LHC) have discovered a new 
gauge boson whose properties are compatible with the 
Standard Model (SM) Higgs boson [1, 2], This discovery 
has initiated a new era of precision studies of the Higgs- 
boson properties. Among the Higgs-boson decay chan¬ 
nels, the decay H —» gg plays an important role in Higgs 
phenomenology. The next-to-leading order (NLO) [3-8], 
the next-to-next-leading order (N 2 LO) [9, 10], and the 
next-to-next-to-next-leading order (N 3 LO) [11] QCD cor¬ 
rections to the total decay width T(H —> gg) under the 
modified minimal subtraction scheme (MS-scheme) have 
been done in the literature. Those improvements provide 
us great chances for achieving precise perturbative QCD 
(pQCD) predictions on T(H —► gg). 

A key problem in making precise pQCD prediction is 
how to set the renormalization scale of the running cou¬ 
pling. Conventionally, the total decay width T(H —> gg) 
is predicted by choosing the Higgs mass Mh as the 
renormalization scale and by varying it within a certain 
range, such as [Mh / 2,2 Mh\, to ascertain its uncertainty. 

However, at any finite order, such a simple choice of 
scale leads to the well-known renormalization scheme- 
and-scale ambiguities. There are also uncanceled large 
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logarithms as well as the “renormalon” terms in the high- 
orders that diverge as (n!/3"a™), which can give sizable 
contributions to the theoretical estimation and dilute the 
pQCD convergence. Schematically, one can express the 
total decay width T(H —> gg) as K x rBorn(H —>■ gg), 
where Tsorn stands for the tree-level or the leading-order 
(LO) decay width. It has been shown that the QCD 
correction factor K under the MS-scheme follows the 
trends [11] 

K m ~ 1 + 0.65 + 0.20 + 0.02 +(1) 

which shows that the NLO and the N 2 LO corrections are 
about 65% and 20% of the Born term, indicating a slow 
pQCD convergence under the MS-scheme. Moreover, as 
will be shown later, there are large scale uncertainties for 
each loop terms under the conventional scale-setting. It 
is thus helpful to find a proper scale-setting approach to 
achieve a more reliable pQCD prediction. 

A guiding principle for resolving the renormaliza¬ 
tion scheme-and-scale ambiguities is that physical results 
must be independent of theoretical conventions. As has 
been argued in Refs. [12-16], if one fixes the renormaliza¬ 
tion scale of pQCD series using Principle of Maximum 
Conformality (PMC), all non-conformal {^}-terms in 
pQCD series are then resummed into the running cou¬ 
pling, and one thus obtains a unique, scale-fixed, and 
scheme-independent prediction at any finite order. The 
PMC has a solid theoretical foundation, satisfying the 
renormalization group invariance [17] and all other self- 
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consistency conditions from renormalization group equa¬ 
tion (RGE) [18]. The PMC provides the general pro¬ 
cedure underlies the commensurate scale relations [19], 
which ensures the scheme independence of the predictions 
under various schemes. A demonstration of the scheme 
independence at any fixed-order has been given with 
the help of Ra-scheme [15, 16], a systematic generaliza¬ 
tion of the minimal subtraction renormalization scheme. 
The PMC reduces in Abelian limit to the standard Gell- 
Mann-Low scale-setting method used in QED [20]. The 
PMC has been successfully applied to many high-energy 
processes, cf. a recent review on the importance of proper 
renormalization scale-setting for QCD testing at high- 
energy colliders [21], which can also be applied to pro¬ 
cesses with multiple physical scales [22, 23]. In the pa¬ 
per, we shall adopt the PMC with the goal of eliminating 
the renormalization scale ambiguity and achieving a more 
accurate pQCD prediction. 

A PMC analysis of the decay H —>• gg under the MS- 
scheme has been done in Ref. [24]. The pQCD conver¬ 
gence of this process cannot be greatly improved as the 
usual PMC applications, especially its N 2 LO contribu¬ 
tion is still about —20% of the total decay width, in¬ 
dicating there is somewhat larger residual scale depen¬ 
dence. However, it is noted that the simple PMC analy¬ 
sis done in Ref. [24] could be inaccurate. This is because 
that the present process involves the three-gluon vertex 
at the lowest order, thus the scale-setting problem should 
be much more involved [25]. Only those {/3j}-terms that 
are pertained to the renormalization of running coupling 
should be absorbed into the running coupling so as to 
achieve the optimal scales at each perturbative order. 
Thus, special treatment, though difficult, should be paid 
for distributing the {/3i}-terms of the process, which is 
lacking in Ref. [24]. 

To avoid such ambiguity of applying the PMC, we shall 
first transform the results from the MS-scheme to the mo¬ 
mentum space subtraction scheme (MOM-scheme) [26- 
29] , which is based on renormalization of the triple-gluon 
vertex at some symmetric off-shell momentum, and then 
apply the PMC scale-setting. Unlike the MS-scheme, the 
physical MOM-scheme carries information of the vertex 
at specific momentum configuration. This external mo¬ 
mentum configuration is non-exceptional and there are 
no infrared issues, thus avoiding the confusion of dis¬ 
tinguishing {/3j}-terms. In the literature, as an alterna¬ 
tion of the MOM-scheme, the minimal momentum space 
subtraction scheme (mMOM-scheme) has also been sug¬ 
gested [30], which is an extension of the MOM-scheme 
on the ghost-gluon vertex and allows the strong running 
coupling to be fixed solely through a determination of 
the gluon and ghost propagators. The mMOM-scheme 
can be related to the MS-scheme at four-loop level, 
thus the four-loop mMOM /3-function can be determined 
with the help of the relations between mMOM and MS- 
couplings [30, 31]. This four-loop mMOM /3-function can 
then be adopted to get the four-loop mMOM-scheme as- 
running behavior and well suit the needs for our present 


four-loop analysis. 

Using the PMC under the MOM-scheme has already 
been suggested by Refs. [32-35] to deal with QCD BFKL 
Pomeron. Those papers show that after applying the 
PMC, the QCD perturbative convergence can be greatly 
improved and the BFKL Pomeron intercept has a weak 
dependence on the virtuality of the reggeized gluon, re¬ 
sulting a much better agreement with the experimental 
data in comparison to the MS-prediction. In those ref¬ 
erences, the NLO transition of the MS-scheme to the 
MOM-scheme has been explained. At present, we shall 
provide a higher-order treatment of such scheme transi¬ 
tion and present the PMC formulas under the mMOM- 
scheme up to four-loop level. 

The remaining parts of the paper are organized as fol¬ 
lows. In Sec. II, we present the technology for calculat¬ 
ing the total decay width T(H —> gg) under the mMOM- 
scheme up to four-loop level. The behavior of the running 
coupling under mMOM-scheme is presented. Numerical 
results and discussions are presented in Sec. III. A brief 
summary will be given in Sec. IV. For convenience, we 
put the PMC coefficients under the mMOM-scheme and 
the Landau gauge in the Appendix. 


II. TOTAL DECAY WIDTH F(H -a gg) UNDER 
THE MINIMAL MOM-SCHEME 

In the literature, the pQCD calculation is usually done 
under the MS-scheme. To apply the PMC, we shall trans¬ 
form the pQCD expressions for the decay H -A gg under 
the MS-scheme into those of the mMOM-scheme by using 
the relation of running coupling between the MS-scheme 
and the mMOM-scheme. 

The MOM-scheme is gauge dependent, and three 
gauges as Landau gauge (£ = 0), Feynman gauge (£ = 1) 
and Fried-Yennie gauge (£ = 3) are adopted in the litera¬ 
ture, where £ stands for the gauge parameter. The ques¬ 
tion is much more involved when the gauge parameter 
£ ^ 0, especially the following suggested extended renor¬ 
malization group method can not be directly applied and 
some alterations must be done to obtain a smooth scheme 
transformation among different running couplings 1 . For 
definiteness, we shall adopt the Landau gauge (£ = 0) to 
do our calculation. 


A. The running coupling and the /3-function to 
four-loop level under the mMOM-scheme 

The scale-running behavior of the running coupling is 
controlled by the renormalization group equation or the 
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where the symbol “A” stands for an arbitrary renormal¬ 
ization scheme, a A = a s , A /4n with a Si A being the strong 
running coupling under the A-scheme. 

By taking the same integral constant as that of 
Refs. [36, 37], the solution of Eq.(2) over the power se- 
(2) ries of 1 /La can be written as, 
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where L A = In (/r 2 /A a ) and A a is the asymptotic scale 
under an arbitrary A-scheme. The {/3^ [[-functions under 
the MS-scheme can be found in Refs. [38-42]. Under the 
Landau gauge, the {/^[[-functions under the mMOM- 
scheme up to four-loop level are [30] 

a™ mom = 11 — 0.66667n/, 

ffMOM = 102 - 12.66667n/, 

ffMOM = 3040.48229- 625.3866771/ + 19.3833n^, 

„9 3 mMOM = 100541.0586 - 24423.3305571/ 

+1625.40224n 2 - 27.49264n^, (4) 

where n/ stands for the number of active flavors. Thus 
the a s -running is fixed if we know the asymptotic pa¬ 
rameter A m MOM well, whose value can be determined by 
using its relation to A^jg. 

It is noted that both 9“ MOM and ^mMOM are the same 
as those of MS-scheme, then, we can apply the extended 
renormalization group method suggested by Ref. [43] to 
evolve the mMOM-running coupling a m MOM(/ir, m MOM) 
at the scale /7 r , m MOM “adiabatically” into the MS- 
running coupling a Mg(My M§) a t sca l e M r MS n0 ^ 0n ^ 
in scale but also in scheme. Thus a more reasonable 
and accurate perturbative expansion of a MS^rMs) over 
a m MOM (+,mMOM) can be achieved [44], As a special 
case, when taking the same arguments for the run¬ 
ning coupling under different schemes, i.e. Mr ms = 
Mr-.mMOM = /+ j we get a perturbative series 

OO 

a Ms(t+) = 1 r '' °mMOM(+)i (5) 

2—1 

whose first three coefficients are 


( 6 ) 

(7) 


2ft ln A mMOM + 4 ^2 ln 2 

Ams A ms 


On the other hand, the coefficients r,; have been di¬ 
rectly calculated up to four-loop level [30, 31, 45-47] 

r 1 = D 1 = 1, (9) 

r 2 = D 2 = -14.0833 + 1.1111171/, (10) 

r 3 = D 3 = -78.7945 + 9.862ti/ + 1.23457n), (11) 

r 4 = D 4 = 862.512 - 328.144 ti/ + 19.256ti^ + 1.37174fcf) 


As a combination of Eqs.(6, 7, 8) and Eq.(9, 10, 11), we 
obtain the wanted relations among the asymptotic scales 
A m MOM and A^jg up to three-loop level, i.e. 


A m MOM 
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(13) 

(14) 

(15) 


where the symbol • • • stands for higher-order equations 
derived through a four-loop and even higher level com¬ 
parisons. The first equation (13) is derived from a 
two-loop comparison, which equals to the one given by 
Refs. [30, 31] that is derived via a different approach. As 
a byproduct, those relations could be treated inversely as 
a consistency check of the complex high-order coefficients 
Di calculated in the literature. For example, the equiva¬ 
lence of Eq.(13) and Eq.(14) can be used as a cross-check 
of the correctness of D 2 and D 3 . 


B. The PMC analysis of H —1 gg under the 
mMOM-scheme up to order a 5 s 

Under the usual MS-scheme, the pQCD prediction for 
the total decay width T(H —> gg) up to n t h-loop level 
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can be written as 


T(H^gg) = J2 C ™ S ^)ai±l(Vr), (16) 
2=1 


where /z r stands for an arbitrary choice of renormaliza¬ 
tion scale, and Cj MS are j t h-loop coefficients under the 
MS-scheme, whose values up to four-loop level can be 
found in Ref. [11]. With the help of Eqs.(5, 9, 10, 11, 12), 
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it can be transformed into the mMOM-scheme as 

n 

T(H -+ gg) = £<«" MOM (Mr) a^ 0 M 

i— 1 

G f M^j i+i / x 

“ 36 ^tt ^ \^o f ) mM ° M 
+0 (a^MOM^r)) ) (17) 

where Gf is the Fermi constant. For convenience, we put 
the mMOM-scheme c+j coefficients with explicit scale de¬ 
pendence in the Appendix. As a further step, following 
the PMC /^-scheme [15, 16], we can transform the n/- 
series into the following {/3™ MOM }-power series, 


T(H -+ gg) = [fi,oaUoM(fr) + (r 2 ,o + 2ft mMOM r 2 ,i)ai M OM(/ i r ) + fa.o + 2/3] nMOM r 2 ,i + 3^ MOM r 3 ,i 

+3(C OM ) 2 ^,2)a 4 mM oM(fr) + (r 4 ,o + 2 /3 2 mMOM r 2jl + 3^ MOM r 3 ,i + 7^ MOM ^ MOM r 3 ,2 
+4/5“ MOM r 4 ,i + 6(/3™ MOM ) 2 r4, 2 + 4(/?™ M °4 3r 4,3)amMOM(Mr-)] + G (a„ MOM (/J r )) . (18) 


The mMOM-scheme conformal coefficients r+o with i = 
(1, 2, 3,4) and the non-conformal coefficients r\ t j with 1 < 
j < i < 4 are also put in the Appendix. 

The non-conformal { / 9‘ nMOM }-terms determine the op¬ 
timal behavior of the running coupling via RGE. Af¬ 
ter applying the PMC, i.e. by absorbing/resumming all 
{/3f lMOM }-terms into the running coupling, the pQCD 
series (18) can be finally simplified into the following 
scheme-independent conformal series: 

C 4 

r(i7 ->■ gg) = X] ri ’° “mMOM(ft) 4 -(19) 

V 2=1 , 


Here Qi stands for the PMC scale at each perturbative 
order, whose values are determined by RGE and are usu¬ 
ally different at different orders. Up to four-loop QCD 
corrections, the LO PMC scale Qi, the NLO PMC scale 
Q 2 and the N 2 LO PMC scale Q 3 are 
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( 21 ) 

( 22 ) 


At the four-loop level, the N 3 LO PMC scale Q 4 is unde¬ 
termined due to unknown higher-order {/3j mMOM }-terms, 
and we set its value as the latest determined PMC scale 
Q 3 . After applying the PMC, there are two kinds of resid¬ 
ual scale dependence. The PMC scales are in perturba¬ 
tive series such that to eliminate all known non-conformal 


/9-terms properly via RGE. The first one is from the PMC 
scales themselves, which are due to “uncalculated” high- 
order /3-terms. The second one is from the final pertur¬ 
bative term, in which we have no /9-terms to set its PMC 
scale and we usually take it as the final PMC scale deter¬ 
mined at one-order lower. Generally, those two residual 












5 


scale dependence are highly suppressed, and the conven¬ 
tional renormalization scale uncertainty to the total and 
individual decay widths at each order can thus be greatly 
suppressed. 


increases with the increment of g, r for /j r ^ 450 GeV, and 
then it slightly decreases with the increment of fi r , which 
changes down to 355 keV for g, r = ITeV. Moreover, for 
the case of g, r = Mjj, Table II shows that the K factor 
under the mMOM-scheme follows the trends 


III. NUMERICAL RESULTS AND 
DISCUSSIONS 



1-loop 

2-loop 

3-loop 

4-loop 

A“ s =5 (GeV) 

0.0904 

0.233 

0.214 

0.214 

A“ m = ° 5 m (GeV) 

0.158 

0.406 

0.373 

0.373 


TABLE I. The asymptotic scales A^ s =5 and A“^T < 5 M at differ¬ 
ent loop-levels under the MS-scheme and the mMOM-scheme. 


To do the numerical calculation, we adopt Gf = 
1.16638 x 10~ 5 GeV“ 2 , the Higgs mass M H = 126 GeV, 
and the top-quark pole mass m t = 173.3 GeV [48]. 
The QCD asymptotic scales are determined by using 
the world average of the running coupling at the scale 
M z , a s (M z )= 0.1185 [49]. The asymptotic QCD scales 
under both the MS-scheme and the mMOM-scheme for 
rif = 5 at different loop-levels are presented in Table I. 
The asymptotic scales with different flavors can be deter¬ 
mined via the usual matching, cf. Ref. [50]. Moreover, to 
be self-consistent, we shall adopt the nth-loop a s -running 
to predict the decay width T(H —> gg) up to nth-loop 
QCD corrections. 


A. Total decay width F(H —> gg) at a 3 -order before 
and after the PMC scale-setting 



LO 

NLO 

N 2 LO 

N 3 LO 

Total 

Ti\M H /2 

284.0 

104.5 

-28.2 

-17.5 

342.8 

L i | 

224.6 

120.4 

15.0 

-9.1 

350.9 

7i\2M H 

184.4 

124.3 

41.9 

5.5 

356.1 

P i\/tM H 

155.4 

123.2 

59.5 

20.1 

358.2 


TABLE II. Total and individual decay widths (in unit: keV) 
of H -a gg at a 3 -order under conventional scale-setting and 
the mMOM-scheme. Ti stands for the individual decay width 
at each order with i = LO, NLO, N 2 LO, and N 3 LO, re¬ 
spectively; r Total = ^Alb stands for total decay width. 
g r = Mh/2, Mh, 2Mh, and 4 Mh, respectively. 

We present the total and individual decay widths of 
H -A gg at a^-order under conventional scale-setting and 
the mMOM-scheme in Table II, in which the results are 
given by varying fi r £ [M H /2,4M H \. r To tai = 
stands for total decay width, where Tj stands for the 
individual decay width at each order with i = LO, NLO, 
N 2 LO, and N 3 LO, respectively. It is found that under 
conventional scale-setting, the total decay width Ti-otai 


-KmMOMlconv. ~ 1 + 0.54 + 0.07 — 0.04. (23) 


In comparison to Eq.(l), it indicates that a more con¬ 
vergent pQCD series than that of MS-scheme can be 
achieved by using the mMOM-scheme even before ap¬ 
plying the PMC scale-setting. After applying the PMC, 
the pQCD series shall be further improved as 

AT m MOM|pMC ~ 1 + 0.30 — 0.10 — 0.04. (24) 


Under conventional scale-setting, a single renormaliza¬ 
tion scale is “guessed”, which, at high orders, may give a 
reasonable prediction for a global observable such as the 
total decay width or total cross section that is close to 
the experimental result; however, the corresponding pre¬ 
dictions for more detailed observables such as the correla¬ 
tions could be inaccurate or even wrong. As an example: 
The two-loop prediction for the total cross-section of the 
top-pair production agrees with the CDF and DO mea¬ 
surements by simply setting the scale as the top-quark 
mass; however one then finds a large discrepancy with 
the top-pair forward-backward asymmetry measured at 
the Tevatron. In contrast, the PMC prediction predicts 
both the total cross-section and the forward-backward 
asymmetry correctly [51-54], which is also confirmed by 
the PMC analysis of top-pair production at the LHC [55] . 

Table II shows that by varying fi r £ [Mjy/2,4M#], 
the total decay width under conventional scale-setting 
shall be changed by about ±2% from its central value at 
fj, r = Mh- This indicates that the scale dependence for 
the total decay width are small at the present af-order. 
By analyzing the pQCD series in detail, it is found that 
the scale errors are rather large for each perturbative 
term P,, thus such a small scale error (±2%) for the total 
decay width is due to cancelations among different orders. 
For definiteness, we define a ratio to show explicitly how 
the individual decay width Tj at each order changes with 
different choices of scale, i.e., 


Hi = 


r *k -r 


i\fi r =M H 




x 100%, 


(25) 


where P^ stands for the individual decay width at each 
order and i = LO, NLO, N 2 LO, and N 3 LO, respectively. 
Under conventional scale-setting, we get 

klo = 26%, knlo = —13%, 

k n 2 lo = —288%, kn 3 lo = 92% (26) 


for n r = Mh/2; we get 

klo = —18%, knlo = 3%, 
k n 2 lo = 179%, kn 3 lo = —160% 


( 27 ) 









6 


LO NLO IsnLO N 3 LO Total 


TilMjj/2 

298.0 

87.9 

-30.0 

-12.5 

343.4 

r i | M h 

298.2 

90.2 

-30.0 

-12.5 

345.9 

r i\iM H 

298.5 

91.9 

-30.0 

-12.5 

347.9 


298.8 

93.3 

-30.0 

-12.5 

349.6 


TABLE III. The PMC predictions for the total and individual 
decay widths (in unit: keV) of H —► gg at o^-order under the 
mMOM-scheme. T; stands for the individual decay width at 
each order with i = LO, NLO, N 2 LO, and N 3 LO, respectively; 
rTotal = stands for total decay width. g, r = Mh/ 2, 

Mir, 2 Mh, and 4 Mh, respectively. 



/ir(GeV) 

FIG. 1. The LO, NLO and N 2 LO PMC scales Qi, Qi and Q 3 
versus the initial choice of scale fj, r for H —> gg at a®- Qrc j er 
under the mMOM-scheme, which are shown by solid, dotted 
and dashed lines, respectively. 

for fj, r = 2 Mh- Those large Ki values indicate that under 
conventional scale-setting, one cannot decide what is the 
exact decay width for each perturbative order. 

As a comparison, we present the PMC prediction for 
the total and individual decay widths under the mMOM- 
scheme in Table III, in which the results are given by tak¬ 
ing /i r = Mh/ 2, Mh 1 2Mh and 4 Mh, respectively. Ta¬ 
ble III shows by varying fj, r G [Mh/ 2,4M#], the total de¬ 
cay width under PMC scale-setting shall only be changed 
by about ±1%. Moreover, the renormalization scale er¬ 
rors for the individual decay widths r, have also been 
greatly suppressed by PMC scale-setting, e.g. except for 
the NLO ratio |«nlo| < 4%, all other Ki are less than 
0.2% by varying fx r within the region of G [M H /2, 4M#]. 
This can be explained by the fact that after applying the 
PMC, we can fix the a s -running behavior, or equivalently 
the renormalization scale, at each perturbative order via 
the using of RGE. 

More specifically, for the present process at the a®- 
order, we have three PMC scales Qi, 2 , 3 , which are de¬ 
termined by Eqs.(20, 21, 22). They are functions of 
the initial scale ^i r . Generally, the initial scale /j, r can 
be chosen arbitrarily, which only needs to be within the 
perturbative region to ensure the reliability of a pQCD 
calculation. The PMC adopts the {/3™ MOM }-terms, de¬ 


termined by RGE, and fixes the optimal scales of the pro¬ 
cess. We present the PMC scales versus the initial scale 
/j, r in Fig.(l). Those optimal PMC scales are smaller 
than the conventional choice of scale (Mh) to a certain 
degree. Fig.(l) shows that the PMC scales Qi, 2,3 are 
highly independent on the choice of /i r , and their values 
are almost fixed to be 

Qi ~ 55GeV, Q 2 ~ 54GeV, Q 3 ~ 109GeV. (28) 

This explains why the individual decay widths at each 
order and the total decay width are highly independent 
to the initial choice of scale. Thus the conventional scale 
ambiguity is cured. 

The {/3™ MOM }-terms are governed by RGE, which are 
generally different at different orders, thus the PMC 
scales Qi at different orders are different from each 
other [13]. The PMC scales themselves are in pertur¬ 
bative series [15, 16], i.e., for H —> gg at order a®, the 
LO PMC scale Q\ is at the accuracy of next-to-next- 
to-leading logarithmic order (NNLLO); the NLO PMC 
scale Qi is at the accuracy of NLLO; and the N 2 LO 
PMC scale Q 3 is at the accuracy of LLO. This way, the 
PMC scales can be improved when more-and-more QCD 
loop corrections are considered, and the unknown higher- 
order {/3f lMOM }-terms shall lead to residual scale depen¬ 
dence. In principal, as shown by Refs. [23, 51-57], due 
to both the a s -power suppression and the exponential- 
suppression, such residual scale dependence is negligible. 
For the present case, we observe that there is a somewhat 
larger residual scale uncertainty in comparison to previ¬ 
ous PMC applications, i.e. a sizable ~ 1 % residual scale 
dependence has been found to total decay width. This 
indicates that there are large contributions from the un¬ 
known {/3™ MOM }-terms. More explicitly, by using the 
formulas ( 20 , 21 , 22 ), we obtain 

In QI/^I^Mh ^ -1.83 + 0.51 + 0.19 + 0(a 3 s ), (29) 
lnQ 2 2 /n% r=MH ~ -2.42 + 2.29 + 0(a 2 s ), (30) 

lnQ 2 3 /nXr=M H =* —0.29 + 0(a s ). (31) 

The PMC scale Q 1 shows a good pQCD convergence as 
the usual PMC applications. However, the pQCD con¬ 
vergence of Qi is questionable, the magnitude of its NLO 
term is at the same order of the LO term, which explains 
the larger residual scale dependence shown by Fig.(l). It 
also explains why there is somewhat larger residual scale 
dependence to Pnlo, as shown by Table III. Thus to fur¬ 
ther improve our PMC predictions, we need to finish even 
higher-order corrections, such as a five-loop calculation, 
so as to achieve an accurate NLO PMC scale Qi and a 
precise prediction of the NLO decay width Tnlo with 
less residual scale dependence. 


B. Total decay width up to different a s -orders 

In above subsection, we have shown the properties of 
total decay width at order a® under the mMOM-scheme. 
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FIG. 2. Total decay width r(H —> gg) versus the scale g r up 
to four-loop level under the mMOM-scheme. The solid, the 
dotted, the dashed and the dash-dot lines are for LO, NLO, 
N 2 LO and N 3 LO total decay width T (H —> gg), respectively. 
The upper diagram is for conventional scale-setting, and the 
lower one is for PMC scale-setting. 


As a step forward, it is helpful to learn how the total 
decay width T(H —> gg) behaves before and after the 
PMC scale-setting when more-and-more loop terms are 
included. In present subsection, we make a detailed dis¬ 
cussion on the properties of total decay width from one- 
loop level to four-loop level. 

We present the total decay width T(H —> gg) un¬ 
der the mMOM-scheme versus the initial scale g r up to 
four-loop level before and after the PMC scale-setting 
in Fig. (2). Fig. (2) shows that under conventional scale¬ 
setting, the LO and NLO total decay width T(H — > gg) 
depend heavily on g r , which then becomes weaker-and- 
weaker when more-and-more loop corrections are taken 
into consideration; and at the four-loop level, the scale 
uncertainty of T(H —> gg) is about ±2% within scale 
region fi r £ [Mh/2,AMh\- This agrees with the con¬ 
ventional wisdom that by finishing a highcr-and-higher 
order calculation, one can get a desirable scale-invariant 
estimate. However, under conventional scale-setting, by 
using a single “guessed” scale, the scale ambiguities and 
scheme-dependence persist at any fixed order: if one uses 
conventional scale-setting for an a™-order pQCD predic¬ 


tion, the large scale ambiguity exists for any known per¬ 
turbative terms, as has been shown by the last subsec¬ 
tion. On the other hand, it is found that after applying 
the PMC scale-setting, the total decay width T(H —> gg) 
with QCD corrections up to LO, NLO, N 2 LO and N 3 LO, 
accordingly, are almost flat versus the scale g r . Fig. (2) 
also shows that after the PMC scale-setting, the value 
of T{H —> gg) shows a faster steady behavior than the 
conventional predictions by including highcr-and-higher 
order corrections, e.g. it quickly approaches its steady 
value with more-and-more loop corrections included. 

As shown by Tables II and III, if setting g r = Mh/2 for 
conventional scale-setting, we get almost the same total 
decay width as the PMC one. This indicates that for con¬ 
ventional scale-setting, the best choice for H —>■ gg should 
be ~ Mh /2 other than the usually suggested Mh- Sim¬ 
ilar choice has been tried for analyzing the gluon-fusion 
channel gg —»■ H [58-60], which has the same topology as 
the decay H gg. In Refs. [58-60], this choice of scale is 
“guessed”, and our present derivation provides a reason 
for this choice 2 . 

It is helpful to estimate the magnitude of “unknown” 
higher-order pQCD prediction. We adopt the way sug¬ 
gested by Ref. [17] for such an estimation, i.e. for a nth- 
loop pQCD prediction 

n 

rr = E C r MOM (^)«jnMOM(^K]), (32) 
2=1 

whose perturbative uncertainty is 

AIT = ± IC MOM (Mr Xmom (QuW\)\max, (33) 

which is calculated by varying g r within the region of 
[Mjj/ 2, 4Mff], and the symbol “MAX” stands for the 
maximum value of |C™ MOM a^ OM | within this scale 
region. Under conventional scale-setting, Qi[g r ] = MU 
under PMC scale-setting, Qi[g r \ are PMC scales and 
QnMOM(^) = GfM/j/( 36V2Tr)r it o are conformal coef¬ 
ficients. This way of estimating “unknown” higher-order 
pQCD prediction is natural for PMC, since after the 
PMC scale-setting, the pQCD convergence is ensured and 
the only uncertainty is from the last term due to unfixed 
PMC scale at this particular order. 

The fixed-order prediction together with its errors on 
the total decay width T(H gg) at different pQCD 
orders before and after the PMC scale-setting are dis¬ 
played in Fig. (3), in which the “Conv.” stands for the 
result under conventional scale-setting. In both cases, 
the predicted error bars from “unknown” higher-order 
terms show a better agreement with their steady points 
when more loop-terms are “known”. 

To compare with the case of conventional scale-setting, 
after the PMC scale-setting, the predicted “unknown” 


2 A detailed PMC analysis to the important Higgs hadro- 

production channel, gg —► H, is in preparation. 
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FIG. 3. Results of total decay width F (H —> gg) = rJ/ ot (in 
unit: keV) at different pQCD orders with their errors AF^° , 
n = (2,3,4), respectively. The crosses and the diamonds 
are for the ones under conventional and PMC scale-settings, 
respectively. The central values are for g r = Mh- 


higher-order contributions are usually smaller than the 
conventional ones [17]. For the present case, there is 
one exception that the predicted AFg ot from a three-loop 
prediction Tg 04 is still large. Numerically it is caused 
by a smaller NLO PMC scale Q 2 (< Mh) and a larger 
N 2 LO conformal coefficient [GfM/j/ (36-^7!") r^o], which 
is about two times larger than the N 2 LO perturbative co¬ 
efficient (7 ™ mom 3 . This exception does not indicate the 
breaking of the general features of the PMC scale-setting, 
which is consistent with the previous observation that the 
NLO PMC scale Q 2 is of large residual uncertainty. In 
fact, for the present decay channel, the predicted error 
bars under conventional scale-setting cannot be treated 
too seriously: the large Ki values, as shown by Eqs.(26, 
27), indicate that one cannot decide the exact values for 
all the individual Ti under conventional scale-setting. 


IV. SUMMARY 

Following the renormalization group invariance, the 
physical observables must be independent on the choice 
of renormalization scheme. The PMC provides a solid 
and unambiguous procedure to set the renormalization 
scale for any QCD process at any finite order. One of 
its key point is to correctly deal with the {/3j}-terms. 
The physical mMOM-scheme is helpful for setting the 
{/?i}-terms of the process unambiguously, especially for 
the processes involving three-gluon or four-gluon vertex 
at lower orders. I11 the paper, we have made a detailed 


3 It is noted that a smaller coefficient (ymMOM j n comparison to 
the conformal one is due to the accidentally large cancelation 
among the conformal and non-conformal terms. 


PMC analysis on the Higgs-boson decay H -A gg within 
the mMOM-scheme up to four-loop QCD corrections. 

It is the first time to apply the PMC to deal with 
N 2 LO and higher order mMOM pQCD predictions. As 
indicated by Table II, the pQCD convergence under the 
mMOM-scheme is better than that of the MS-scheme 
even before the PMC scale-setting. After applying the 
PMC, the pQCD convergence can be further improved. 

Table II shows that the scale dependence for total de¬ 
cay width are small at the present a^-order even un¬ 
der conventional scale-setting, i.e. by varying £ 
[Mh/2,AMh], the total decay width shall only be 
changed by ~ ±2%. Under the conventional scale¬ 
setting, by taking Mh = 126 GeV, we obtain 

r(if —> 99 ) |mMOM,Conv. = = 350.9tg f keV, 

(34) 

where the central value is for = Mh and the er¬ 
ror is for fi r £ [Mh/2,4Mh\- As we have pointed 
out in the body of the text, the scale errors for each 
loop terms are quite large, and a small net scale er¬ 
ror for the four-loop total decay width is due to the 
large cancelations among different orders. More explic¬ 
itly, the scale errors for the decay widths Tj can be 
reexpressed by the parameter Ki (defined in Eq.(25)), 
which are klo £ [26%, — 31%], knlo £ [—13%,2%], 
kmto £ [-288%, 297%] and k N 3 LO £ [92%, -321%] for 
Hr £ [M h / 2 , 4 Mh]- 

After applying the PMC scale-setting, the scale errors 
for either the total decay width T(H -A gg) or its indi¬ 
vidual contributions T, are largely suppressed. By taking 
Mh = 126 GeV, we get 


T(H -a 5 g)| mM OM,PMC = 345.9ti;£ keV, (35) 

where the central value is for fi r = Mh and the error is 
caused by varying \i r £ [Mh/2,AMh\- The central decay 
width is lowered by about 1% in comparison to conven¬ 
tional one. The scale error for the decay width Tj of each 
loop is greatly suppressed, e.g. except for the NLO ra¬ 
tio |knlo| ~ 4%, all other ratios K, are less than 0.2% 
by varying fj, r £ [Mh/2,4Mh]- It should be pointed out 
that the somewhat larger residual scale error as shown by 
Eq.(35) in comparison to previous PMC examples are due 
to the uncalculated high-order /3-terms. As shown by Ta¬ 
ble III, such large error is from the residual scale depen¬ 
dence of the NLO decay width TnlOj whose PMC scale 
Q 2 has a poor pQCD convergence as shown by Eq.(30). 
Thus to further improve the accuracy of the PMC pre¬ 
dictions, one needs to finish even high-order corrections 
such as a five-loop calculation to the present channel so 
as to achieve an accurate Q 2 and a precise prediction of 
Tnlo with less residual scale dependence. 

Finally, if setting the Higgs mass as the recently mea¬ 
sured one by the ATLAS and CMS collaborations, the 
PMC prediction of the total decay width is 

T(H -A gg) Umom.pmc = 339.3 ± 1.7keV, (36) 
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where the first error is caused by taking the Higgs 
mass M h = 125.09 ± 0.21 ± 0.11 GeV [61], and the 
second error is caused by varying /z r £ [M#/2, AMh\- 
To compare with Eq.(35), it shows that a change of 
Higgs mass by 1 GeV, the total decay width shall be 
changed by about 6 keV. Thus the eliminating of the 
renormalization scale error shall inversely help us to get 
an accurate prediction on Higgs mass, and etc.. 
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APPENDIX: THE COEFFICIENTS aj AND nj 
UNDER THE MMOM-SCHEME 


As mentioned in the body of the text, the Cjj coeffi¬ 
cients under the mMOM-schemes can be obtained from 
those of MS-ones with the help of Eq.(5). The mMOM 
Cij coefficients under Landau gauge with full renormal¬ 
ization scale dependence are 


ci,o = 16, 

c 2 ,0 = 1069.33- 352 In 


m 2 h 


M 


c 2 i =-39.1111+ 21.3333 In 

M 2 

M 2 M 2 

c 3 ,o = 31202.1 + 608 In —if - 38552 In —f 


(37) 

(38) 

(39) 


+5808 In 


2 M 2 h 


AT 


(40) 


M 2 M 2 

c 3 ! = -4043.54 + 170.667 In —£ + 3834.67 In —f 

mi at 


Mr 


(41) 


Ml 2 Ml, , 

c 3 2 = 41.2236 - 78.2222 In —f + 21.3333 In 2 -^-(42) 

Mr Mr 

M 2 M 2 

c 4 ,o = -88214.8 + 6688 In 2 —f + In —f-(34008.9 


o 

mi 


mi 


-26752 In ^%) - 1.7974 x 10 6 In 

M r Mr 

+902000 In 2 - 85184In 3 ^K, (43) 

Mr Mr 

M 2 M 2 

c 4>1 = -128899 + 1472 In 2 —f + In —f-(8037.93 


Ml x M 2 

-5888 In -ML) + 333737 In —f 


m2 


“m 


3 1V *H 


2 Mrr 

-145717 In 2 -ML + 15488 In' „ 
Mr Mr 


(44) 


2 Ml Ml, 

c 4 2 = 3663.24 - 113.778 In 2 —f + In —^(53.3333 

M 2 M 2 

+455.111 In ^f)- 14703.1 In MJL 
Mr Mr 


+7239.11 In 2 —if. — 938.667 In 3 (45) 


Mr 


m: 


Ml 


c 4 3 = 9.97767 + 109.929 In -ML - 104.296 hT -ML 

Mr Mr 

, Ml 

+18.963 In 3 — LL' (45) 

Mr 

There are two types of logarithmic terms In Ml/[H r 
and In m 2 //x 2 in those expressions, and we have used 
In Ml/mi to replace (In Ml/gi — In mi / ni) for brevity. 

The expansion of the mMOM-scheme coefficients r^j 
up to four-loop level under the Landau gauge are 


1+0 = Ci, 0 , 

3C2 1 

^ = 2^T’ 

33c2,i 

1-2,0 = C2,0 + 


2 ’ 


r 3,2 = 


9c 3 ,2 


1 - 3,1 = - 

1 


2n(l + n) ’ 

6(c 3 ,i + 33 c 3j2 ) - H 4 c 2 ,i 


(47) 

(48) 

(49) 

(50) 

( 51 ) 


4(1 + n) 

1-3,0 = q (8c 3i o + 66(2c 3 ,i + 33c 3 ,2) — 1284 c2,i), (52) 

81c4, 3 


1 - 4,3 = -- 


(53) 


4 n (2 + 3 n + n 2 ) ’ 

1+2 = 0/1 1 \ 2 (o, . ( 989 c 2 i(l + n) + 18((2c 4 ,2 

8(1 + n) z ( 2 + n) 

+99c 4 j3 )(1 + n) — 38 c 3j2 (3 + 2 n)) 
+48c 2 ,i(1+h)C3), (54) 

i* 4 , 1 = q9/1 , w, , (~14554 c 2 ,i(1 +n) 

32(1 -|- TljyZ Tl) 

—6(2(4 c4 ? i + 33 ( 4 c 4 7 2 99c4 ?3 ))(l + u ) 

- 152 c 3 ,i (1 + n ) + c 3?2 (—8868 - 7584 n)) 

—3162c 2 ,i(1 + h)C 3 ), (55) 

1-4 0 = ttt( 3044 c 2 1 + 8 ( 8 c 4 0 + 132 c 4 1 + 2178 c 4 2 
64 

+35937c 4 , 3 - 1284c 3j i - 42372c 3>2 )), (56) 

where £ 3 is Riemann zeta-function, and n stands for the 
LO a s -power, which equals to 2 for the decay H -+ gg. 
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